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Abstract. 

The orbital degrees of freedom are of vital importance in explanation of various phenomena. 
Among them is the orbital-selective Mott transition (OSMT), which is thought to occur in 
several materials as Ca2-a : Sr a: Ru04 and La„+iNi n 03 n +i. OSMT is usually studied in the 
infinite-dimension limit, and for the time being, it is not clear if it could survive in one- 
dimensional (ID) case. There exist two scenarios for the OSMT: upon increasing the interaction 
in a two-band system i) one of the bands becomes insulating, while the other remains metallic 
and ii) one of the bands becomes empty, while the other may eventually undergo a single- 
band Mott insulator transition. In this work, we present a preliminary study of the two-orbital 
Hubbard model by means of Density Matrix Renormalization Group in ID at quarter-filling, 
where the second scenario seems to be realized. In particular, we study the orbital densities, 
double occupancies and form-factors also in the case of finite inter-orbital inter-site hopping. 



1. Introduction 

The crucial role of the orbital degrees of freedom in explaining several intriguing phenomena 
such as metal-insulator transition and magnetism in doped phthalocyanines as well as the 
superconductivity in heavy-fermion compounds has been recently recognized. On the theoretical 
side, the study of multi-orbital systems is extremely complicated and has been mainly 
accomplished in the infinite-dimension limit by means of the Dynamical Mean-Field Theory, 
although, in these systems, the spatial correlations are essential. 

Another fascinating phenomenon inherent to the multi-orbital systems is the orbital selective 
Mott transition (OSMT). In such a transition, by tuning the strength of the Coulomb repulsion, 
it is possible to open a gap in some band(s), while leaving the other (s) gapless. It is thought 
that OSMT occurs in materials as Ca2_ a: Sr. r Ru04 [HE] and La n+ iNi n 03 n+ i [31 HIE]. 

The simplest multi-orbital model is the two-orbital Hubbard model, whose most general 
Hamiltonian reads as follows: 





Figure 1. Left panel: [/'-dependence of the average occupation number per orbital per site n\ 
and 77-2 at different values of inter-orbital hopping t 12 . U = 10, t 11 /t 22 = 0.5, DMRG 50 sites. 
Right panel: The same for the average double occupancy D\ and D2. 
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Here t a @ parametrizes intra- (t 11 and t 22 ) and inter- orbital (t 12 ) hoppings between nearest- 
neighbor (NN) sites, U and U' are the on-site inter- and intra-orbital Coulomb repulsions, 
respectively, while J is the direct exchange interaction. Rotational symmetry imposes that 
U = U' + 2J. Finally, cj creates an electron with spin a in the a-orbital of site i, 

^ia = \ Yl a ab c laa c iab ^ s ^ ne electron spin operator in the a-orbital of site i, being o k ah the 
Pauli matrices, and the pair operator pi a and the electron density operator rii aa are defined as 

follows: Uiaa = 4aa C iaa' Pia = Ci a ^C ia ^. 

In the present manuscript, we use finite-size Density Matrix Renormalization Group (DMRG) 
method [6] to investigate the zero-temperature properties of the Hamiltonian ([2]). We consider 
chains of 50 sites with open boundary conditions. In the DMRG decimation procedure, we retain 
up to M = 300 lowest eigenstates of the density matrix, which amounts to have a truncation 
error on the sum of the density matrix eigenstates not larger that 10~ 5 . 

2. Results 

It is well known that the (OSMT) occurs at commensurate fillings (n = 1, n = 2 electrons 
per site). In particular, in one-dimensional (ID) case, by means of the product-wave-function 
renormalization-group technique, the existence of two Mott states in the Hamiltonian ([2]) at 
quarter- and half- filling for t 12 = 0, U = 10 and t 11 /t 22 = 0.5 has been shown in Ref. [7]. It was 
found therein that the quarter-filling case corresponds to a complete narrow orbital suppression 
(NOS), while in the half-filling case there appears a true OSMT state. In the present manuscript 
we concentrate on the case of a quarter-filled system (n = 1) and explore the evolution of NOS 
upon introducing a finite inter-orbital hopping t 12 . We have also studied the t 12 = case, already 
analyzed in Ref. [7], in order to compare our finite-size algorithm with their bulk one, as well 
as to have our own t 12 = reference line. In Fig. [H we present our results for average electron 




Figure 2. Fourier transforms of the density (Left panel) and spin (Right panel) static form- 
factors at initial (V = 0) and final (V = 20) points from Fig. [1] and different values of inter- 
orbital hopping t 12 : t 12 = 0, 0.25, 0.5. U = 10, t u /t 22 = 0.5, DMRG 50 sites. 



densities and double occupancies per orbital as functions of U' at U = 10 and J = (U — U')/2. 
For t 12 = 0, NOS occurs approximately at U' c = 9.4 in qualitative agreement with Ref. (7j, 
although in the region of U' > U' c ( J < 0) there is a slight re-population of the otherwise empty 
narrow orbital. Such a re-population might be explained as the effect of the negative J Hund's 
exchange: the system gains more energy by populating both orbitals rather than emptying 
one of them. We have checked, by increasing M (not shown), that this re-population is not a 
DMRG artefact, and the small discrepancy between our results and those of Ref. [7] remains to 
be investigated. Our analysis confirm the conclusion of Ref. [7] that the quarter-filled transition 
is due to NOS. NOS influences inevitably the behavior of the average double occupancy per 
site -Dj (i = 1, 2), as shown in the right panel of Fig[TJ As U' grows for U' < U' c , D\ decreases 
monotonically due to the depletion of the narrow orbital. At the same time, Z?2 grows due 
to the opposite reason and reaches the value D2 ~ 0.02455 at U' c which is very close to the 
Bethe-ansatz value {D BA = 0.025 at half-filling and U = 10, Ref.0). 

In order to have a deeper insight into the nature of the initial and final states involved in 
this transition, we have measured the static density and spin form-factors. As shown in FigJ21 
in the initial state (V = 0), there is a strong density correlation peaked at the wave- vector it, 
meaning that the density has staggered correlations. These staggered correlations are gradually 
suppressed at V > U' c . In such a strong-coupling regime, the wide band becomes uniformly 
occupied leading to an almost featureless density form-factor. The opposite picture is observed 
in the spin channel. At U' = and t 12 = 0, the spin form-factor exhibits a peak at 7r/2, implying 
the enhanced staggered correlations with the period of four lattice spacings, while at t 12 7^ 
the spin correlations are heavily suppressed. In the strong-coupling regime at U' > U' c , only the 
wide orbital is active and the system resembles very much a single-band Hubbard model at half- 
filling. In such a case, large on-site repulsion induces strong antiferromagnetic NN correlations, 
as can be seen from Figj2j 

Upon switching on a finite inter-orbital hopping t 12 , the NOS becomes incomplete, as shown 
in the left panel of Fig. [TJ This can be understood by noting that the electrons in the narrow 
orbital can now reach the NN site also by making the virtual hopping through the wide orbital 
according to the following path: (i, 1) -4 (i + 1, 2) — > (i, 2) — > (i + 1, 1), where (i, a) stands for 
a-orbital on site i. Such an extra possibility effectively widens the narrow orbital. To better 
understand this widening, we have investigated the NOS when t 12 = and at different values of 
t , as shown in the left panel of Fig. [3j We have observed the full NOS for t 11 < 0.6. This full 
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Figure 3. Left panel: [/'-dependence of the narrow band occupation at t 12 = and different 
values of t . In the inset, U' c is shown as function of t . Right panel: The narrow band 
occupation as function of t 12 as the latter goes to zero at several values of U' and t . 

NOS occurs at some U' c (t n ) (shown in the inset) above which n\ starts growing nearly linear. 
This situation changes drastically when t 12 ^ 0. The NOS in this case is never complete even 
for small t 12 , as shown in Fig. [TJ In the right panel of Fig. Owe present a study of the narrow 
band occupation as t 12 — > for different values of U' and t . We have chosen these values in 
such a way that there is a complete NOS at t 12 = 0. It can be seen from Fig. [3] that as soon as 
t 12 is finite, ni is never suppressed completely. 

Summarizing, by using finite-system DMRG, we have studied the influence of the inter-orbital 
hopping t 12 on NOS in quarter-filled two-orbital Hubbard model ([2]). We have found that t 12 
prevents the narrow orbital from complete suppression, by providing an additional contribution 
to the intra-orbital hopping. We have found that at t 12 = 0, a complete NOS occurs at U' c (t 11 ) 
for t 11 < 0.6 while for greater values of t 11 NOS is incomplete for all values of U' € [0,20]. 
On the contrary, our data suggest that a finite t 12 prevents complete NOS even in case of very 
narrow band (t 11 = 0.1). In the regime with U' < U' c we find a strong enhancement of staggered 
charge correlations, while in the case when U' > U' c the staggered spin correlations are favoured. 
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